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Closed-loop  systems  have  been  developed  for  controlling  the  flow  above  a  three-dimensional 
turret.  The  top  of  the  turret  is  hemispherical,  houses  a  flat  optical  aperture,  and  can  rotate 
about  two  axes  (pitch  and  yaw).  The  extent  of  separation  and  concomitant  turbulence  levels 
in  the  flow  above  the  aperture  change  as  the  turret  rotates.  Suction  jet  slots  circumscribing 
the  aperture  serve  as  control  input;  an  array  of  pressure  sensors  on  the  turret  surface 
provides  the  controller  with  information  about  the  state  of  the  flow  above  the  surface.  The 
control  objective  is  to  minimize  the  separation  and  turbulence  in  the  dynamic  environment 
created  by  the  articulating  turret.  The  closed-loop  control  systems  include  dynamical  and 
measurement-based  estimators,  regulators,  Alters,  and  compensators.  These  components 
are  developed  using  both  computational  and  experimental  data,  and  the  control  systems 
are  evaluated  through  a  series  of  control-in-the-loop  CFD  simulations  and  wind  tunnel 
runs.  Controller  designs  and  computational  tests  are  described  in  “Feedback  Flow  Control 
for  a  Pitching  Turret  (Part  I),  and  the  follow-on  wind  tunnel  tests  are  described  here. 

I.  Introduction 

Turbulence,  random  density  fluctuations,  and  aero-optical  distortion  in  the  flow  around  a  turret  can 
degrade  the  performance  of  a  laser  system  housed  within  the  turret.  Laser  light  passing  through  the 
turbulent  region  is  refracted  by  variations  in  density  and  optical  path  differences,  scattering  the  light  and 
reducing  the  laser  intensity.  One  solution  to  the  problem  is  to  adapt  the  optical  system  in  a  way  that 
compensates  for  the  flow  distortion — without  altering  the  flow  itself;  for  example,  a  morphing  lens  that  keeps 
the  laser  focused  in  the  presence  of  aero-optic  distortions  through  contortions  of  the  deformable  mirror.^ 
Another  solution  is  to  control  (minimize) — either  actively  or  passively — the  flow  fluctuations  that  cause 
aero-optical  distortion. 

In  a  passive  flow  control  experiment,  vortex  generators  mounted  on  the  upstream  side  of  a  two-dimensional 
turret  successfully  reduced  distortions  at  certain  aperture  angles  of  attack.^  In  another  experiment,  a 
spanwise  array  of  synthetic  jet  actuators  actively  reduced  flow  fluctuations  and  delayed  separation  on  a 
three-dimensional  turret.^  In  an  open-loop  control  experiment,  synthetic  jets  were  used  to  reduce  the  root- 
mean-square  (rms)  values  of  fluctuating  velocity  in  the  turbulent  flow  around  a  three-dimensional  turret.^ 
Closed-loop  control  systems  are  designed  to  improve  performance  over  open-loop  systems  through  ro¬ 
bustness,  the  ability  to  compensate  for  uncertainties  over  a  range  of  flight  conditions,  and  by  achieving  more 
with  less — reducing  actuator  power,  size,  and  weight  requirements.  A  closed-loop  system  utilized  synthetic 
jets  to  delay  the  onset  of  stall  on  a  NACA  4412  airfoil,  consuming  less  power  than  an  open-loop  system  that 
increased  stall  angle  by  the  same  amount.^’®  More  recently,  a  proportional  feedback  controller  was  used  to 
modulate  synthetic  jets  on  a  3D  turret  and  reduce  fluctuation  levels  in  the  wake.^  The  tests  were  performed 
at  Mach  0.3  in  the  Subsonic  Aerodynamic  Research  Laboratory  (SARL)  wind  tunnel  at  Wright-Patterson 
Air  Force  Base,  OH. 
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^Research  Scientist,  Clear  Science  Corp.,  663  Owego  Hill  Road,  Harford,  NY  13784-0233. 
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Here,  closed- loop  flow  control  tests  have  been  performed  with  a  3D  turret  at  Mach  0.1.  Compressible 
effects  are  not  present  at  this  speed,  but  systems  for  controlling  velocity  fluctuations  in  incompressible 
flows  will  be  effective  in  controlling  density  fluctuations  at  higher  speeds  because  of  the  strong  correlation 
between  velocity  and  density  fluctuations.®  Both  computational  and  experimental  tests  were  performed.  The 
controller  designs  and  computational  tests  are  described  in  “Feedback  Flow  Control  for  a  Pitching  Turret 
(Part  I).”  The  follow-on  experimental  tests  are  described  here. 

II.  Problem  Description 

The  control  problem  involves  fully  turbulent,  incompressible  flow  past  a  three-dimensional  (3D)  turret, 
and  the  control  objective  is  to  minimize  levels  of  separation  and  velocity  fluctuations  above  the  aperture. 
Free  stream  conditions  in  the  computational  simulations  and  the  wind  tunnel  conditions  in  the  experiment 
are  approximately  equal.  The  nominal  free  stream  Reynolds  number  based  on  the  diameter  of  the  turret 
is  450,000,  and  the  free  stream  Mach  number  is  0.1.  The  computational  model  of  the  turret  and  the  wind 
tunnel  test  article  also  match — with  some  differences  in  the  details.  The  turret  consists  of  a  six-inch-diameter 
hemisphere  mounted  on  a  six-inch-diameter  cylinder.  The  cylinder  is  four  inches  in  height.  A  flat  aperture, 
2.8  inches  in  diameter,  is  located  at  the  top  of  the  turret.  The  top,  hemispherical  section  can  rotate  within 
the  cylindrical  stand  in  two  directions:  pitch  and  yaw.  Rotation  about  the  pitch  axis  alone  is  examined  here. 

The  pitch  angle  is  defined  as  the  angle  between  a  vector  normal  to  the  aperture  and  the  free  stream 
vector.  When  the  normal  is  aligned  with  the  turret  centerline,  the  pitch  angle  is  90°.  As  the  flat  aperture 
pitches  back,  it  becomes  a  backward  facing  ramp,  separation  levels  increase  as  the  pitch  angle  increases,  and 
this  increases  levels  of  velocity  fluctuations  in  the  flow  above  and  downstream  of  the  aperture.  The  objective 
is  to  manage  and  minimize  the  separation  and  concomitant  fluctuations  in  the  aperture  field  of  view  as  the 
pitch  angle  is  varied  sinusoidally. 

Control-in-the-loop  CFD  simulations  were  used  to  develop  and  evaluate  control  designs  and  methods 
that  are  then  applied  in  the  wind  tunnel.  The  CFD  and  wind  tunnel  controllers  are  comprised  of  the 
same  components  (described  in  “Feedback  Flow  Control  for  a  Pitching  Turret  (Part  I)”);  however,  these 
components  are  calibrated  with  CFD  data  for  the  control  simulations  and  are  calibrated  with  experimental 
data  for  the  control  tests  in  the  wind  tunnel.  The  intent  is  not  to  use  experimental  data  to  validate  CFD 
models  and  methods  but  rather  to  develop  and  evaluate  closed-loop  control  designs  for  the  wind  tunnel 
through  CFD  simulations. 

III.  The  Experimental  Test  Configuration 

The  test  facility  at  Syracuse  University  is  a  Gottingen-type,  closed,  recirculating  tunnel  in  a  horizontal 
configuration  with  continuously  variable  speeds  of  less  than  4  m/s  to  approximately  70  m/s.  The  test  section 
is  constructed  of  optical  Plexiglas  with  dimensions  of  61  cm  (width)  x  61  cm  (height)  x  2.4  m  (length).  The 
turret  is  mounted  directly  to  the  floor  of  the  tunnel.  The  left  panel  of  Figure  1  contains  a  photograph  of 
the  test  article,  and  the  right  panel  contains  a  schematic  of  the  top  of  the  test  article.  Two  concentric  rings 
of  suction  jet  slots  circumscribe  the  aperture.  Slots  on  one  side  of  the  centerplane  are  tied  to  a  common 
manifold,  and  slots  on  the  other  to  another  manifold.  Each  manifold  is  connected  to  a  vacuum  system  with  a 
valve  located  between  it  and  the  slots.  A  total  of  30  pressure  sensors  are  located  on  and  around  the  aperture. 
Surface  pressures  are  sampled  at  a  rate  of  10  kHz  using  a  NATIONAL  INSTRUMENTS  PXI-based  800  MHz 
signal  conditioner  with  24-bit,  high-resolution  A/D  converters  and  anti-aliasing  filter.. 

Suction  jet  velocity  is  controlled  by  cycling  the  suction  valves  at  25  Hertz  and  modulating  the  duty  cycle. 
Velocity  measurements  were  taken  just  above  the  slots  using  a  hot  wire  probe  to  determine  the  relationship 
between  suction  velocity  and  duty  cycle.  Erom  the  left  panel  of  Eigure  2,  the  relationship  is  approximately 
linear  between  duty  cycles  of  30%  and  90%.  The  right  panel  of  Figure  2  contains  a  schematic  of  the  Particle 
Imaging  Velocimetry  (PIV)  system  that  is  used  to  measure  velocity  in  the  flow.  The  DANTEC  PIV  System 
measures  two  components  of  velocity  in  the  centerplane  at  a  rate  of  4  kHz  using  a  DANTEC  FLOWMAP 
system.  The  velocity  measurements,  surface  pressure  measurements,  and  suction  valve  duty  cycle  are  phase 
aligned  or  synchronized  in  time. 

Figure  3  contains  photos  of  the  mechanism  that  effects  rotation  about  the  pitch  and  yaw  axes,  housed 
under  the  top  of  the  turret.  A  drive  shaft  connects  a  stepper  motor  to  the  turret  top  through  a  gearing 
system,  generating  rotation  about  the  pitch  axis.  Another  stepper  motor  generates  yaw  rotation.  The  motors 
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Figure  1.  The  test  article  (left)  and  a  schematic  of  the  top  of  the  turret  with  sensor  and  actuator  arrays(right). 
Jet  slots  are  marked  as  lines  and  sensors  as  open  circles. 


Figure  2.  Suction  jet  velocity  versus  suction  valve  duty  cycle  (left)  and  a  schematic  of  the  PIV  window  (right). 


are  rated  to  58.61  Newton-cm  of  torque  and  are  controlled  by  the  data  acquisition  system,  Lab  VIEW. 

IV.  Static  and  Dynamic  Pitching 

The  first  series  of  wind  tunnel  tests  examine  pitch  angles  between  110°  and  120°.  These  tests  were 
conducted  with  no  actuation  and  with  constant  actuation  at  30%,  50%,  70%,  and  90%  duty  cycles  (open 
loop).  The  turret  is  stationary  at  each  individual  pitch  angle  in  a  set  of  static  runs,  and  the  turret  rotates 
over  the  pitch  angle  range  in  dynamic  tests.  The  intent  is  to  examine  the  effect  of  the  turret  motion  on 
the  flow.  Velocity  is  sampled  at  4  Hertz  and  surface  pressures  at  10  kilo-Hertz.  The  top  panels  of  Figure  4 
contain  contours  and  streamlines  of  time-averaged  velocity  in  the  PIV  centerplane  at  static  pitch  angles  of 
115°  and  120°  (with  no  actuation).  In  both  cases,  the  flow  is — on  average — fully  separated  in  the  region 
above  the  aperture.  The  pitch  angle  of  115°  marks  the  point  of  incipient  separation  with  a  stationary  turret: 
flow  was  separated  above  the  aperture  in  one  run  at  this  angle  and  attached  in  another.  The  flow  was  always 
attached  at  lower  static  angles  and  always  separated  at  higher  angles.  The  bottom  panels  of  Figure  4  contain 
time  histories  of  surface  pressure  from  five  sensors  located  on  the  turret  centerline:  Sensor  Nos.  1,  14,  21, 
28,  and  7  (see  the  sensor  locations  in  the  right  panel  of  Figure  1). 
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Figure  3.  Photos  of  the  two-degree-of-freedom  rotational  mechanism  housed  inside  the  turret. 
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Figure  4.  Contours  and  streamlines  of  time-averaged  velocity  on  the  PIV  centerplane  (top  panels)  and  time 
histories  of  surface  pressure  from  five  sensors  located  on  the  turret  centerline  (bottom  panels)  at  static  pitch 
angles  of  115°  (left)  and  120°  (right)  with  no  suction  jet  actuation. 


Next,  the  turret  pitches  dynamically  over  a  range  from  110°  to  120°  using  the  following  prescribed 
function: 

6»(t)  =  [115  -  5cos(wi)]  yfg  (1) 

where  co  =  dO/dt  =  67r/25  rad/s.  Therefore,  one  pitching  period  is  8g  seconds.  Again,  velocity  is  sampled 
at  4  Hertz  and  pressure  at  10  kilo-Hertz. 
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Figure  5  contains  results  from  the  dynamic  run  with  no  actuation.  The  data  has  been  ensemble  averaged 
over  a  number  of  pitching  cycles.  Comparing  the  top  and  bottom  panels  on  the  left,  the  flow  is  attached  as 
the  pitch  angle  increases  through  115°  and  detached  as  it  decreases  through  115°.  In  fact,  when  the  angle 
is  increasing,  the  flow  remains  attached  until  119° — at  which  point,  it  becomes  massively  separated.  As 
the  angle  reaches  the  maximum  of  120°  and  begins  to  decrease,  the  flow  remains  separated  until  112°  when 
it  re-attaches.  The  bottom  right  panel  of  Figure  5  contains  time  histories  of  surface  pressure  from  Sensor 
Nos.  1,  14,  21,  28,  and  7.  The  fluctuating  pressure  ahead  of  the  aperture  (Sensor  No.  1  in  the  topmost 
time  series)  increases  as  turret  pitches  back  and  the  flow  separates.  However,  the  sensors  located  on  and 
downstream  of  the  aperture  (Sensor  Nos.  14,  21,  28,  and  7  in  the  lower  four  time  series)  exhibit  an  opposite 
trend:  fluctuating  pressures  decrease  as  the  flow  separates.  Comparing  results  from  the  static  and  dynamic 
runs  reveals  that  the  influence  of  dynamic  rotation  is  significant  even  when  the  pitching  rate  is  modest. 
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Figure  5.  Velocity  contours  and  streamlines  in  the  PIV  centerplane  from  the  run  with  dynamic  pitching  and 
no  actuation  at  dynamic  pitch  angles  of  115°  (increasing)  in  the  top  left  panel,  120°  in  the  top  right  panel,  and 
115°  (decreasing)  in  the  bottom  left  panel.  Time  histories  of  five  surface  pressures  on  the  turret  centerline 
from  the  dynamic  run  in  the  bottom  right  panel. 


Open-loop  flow  control  runs  were  executed  with  a  pitching  turret  and  steady  suction  at  30%,  50%,  70%, 
and  90%  duty  cycles.  Figure  6  contains  results  from  the  dynamic  run  at  50%  duty  cycle.  Again,  the  data 
are  ensemble  averaged  over  a  number  of  pitching  cycles.  From  the  panels  containing  velocity  contours  and 
streamlines,  steady  suction  at  the  50%-duty-cycle  level  moves  the  separated  flow  downstream  and  below 
the  aperture.  Comparing  the  bottom  right  panels  of  Figures  5  and  6,  steady  suction — and  the  elimination 
of  separation  above  the  aperture  that  results — actually  increases  pressure  fluctuations  in  all  but  one  of  the 
sensors.  Assuming  that  separation  promotes  turbulent  fluctuations  in  the  flow,  the  only  sensor  that  follows 
the  expected  trend  is  Sensor  No.  1.  It  is  also  the  only  sensor  located  upstream  of  the  jet  slots.  Root-mean- 
square  values  of  fluctuating  surface  pressures  are  higher  than  baseline  values  in  all  of  the  open-loop  runs  (for 
all  of  the  sensors  but  Sensor  No.  1).  The  same  trend  was  confirmed  at  a  duty  cycle  of  100%  (valves  fully 
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open) — eliminating  transmission  of  vibrations  from  the  valves  as  a  possible  source  of  the  trend. 
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Figure  6.  Velocity  contours  and  streamlines  in  the  PIV  centerplane  from  the  run  with  dynamic  pitching  and 
open-loop  control  (50%  duty  cycle)  at  dynamic  pitch  angles  of  115°  (increasing)  in  the  top  left  panel,  120°  in 
the  top  right  panel,  and  115°  (decreasing)  in  the  bottom  left  panel.  Time  histories  of  five  surface  pressures  on 
the  turret  centerline  from  the  dynamic  run  in  the  bottom  right  panel. 


V.  Initial  Closed-loop  Control  Run 

The  initial  closed- loop  control  run  in  the  wind  tunnel  utilized  a  single  pressure  sensor  (Sensor  No.  1), 
chosen  because  its  rms  values  decrease  as  duty  cycle  increases  in  the  open-loop  tests.  A  simple  proportional 
controller  was  employed  in  the  initial  test.  The  control  objective  is  to  minimize  fluctuations  in  the  single¬ 
sensor  pressure  measurement,  using  suction.  The  duty  cycle  of  the  suction  valve  is  modulated  by  the 
controller.  The  turret  pitches  at  a  maximum  rate  of  3.768  degrees  per  second  between  110°  and  120° 

according  to  the  function  in  (1).  One  pitching  period  is  8^  seconds.  Velocity  is  sampled  at  4  Hertz  for  a 

total  of  102  pitching  cycles  (850  seconds),  yielding  3400  PIV  snapshots.  This  corresponds  to  100  snapshots 
over  every  three  pitching  cycles,  and  PIV  data  are  distilled  by  averaging  values  from  34  three-cycle  periods. 
Pressure  is  sampled  at  10  kilo-Hertz.  All  of  the  20  jet  slots  were  activated,  operating  as  a  single  unit. 

Four  control  runs  were  executed,  using  four  different  types  of  pressure  filters: 

•  Band  Pass:  Frequencies  (Hz)  in  the  range  [30,  500] 

•  Band  Stop:  Frequencies  (Hz)  in  the  range  [0,8]  U  [12,500]. 

•  Low  Pass:  Frequencies  (Hz)  in  the  range  [0,500]. 

•  Very  Low  Pass:  Frequencies  (Hz)  in  the  range  [0,8]. 
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Performance  output  is  defined  as 


^  /  Npl^/  , 

-  y^A[uiix^,t)  -Uiix,)]"^  +  [u2iXi,t)  -U2iXi)]‘^} 


^  tVpiv 


\  i=i 


(2) 


where  ui,2(£Ci,t)  are  the  instantaneous  velocities  in  the  centerplane,  and  ui^2{xi)  are  the  mean  (time- 
averaged)  velocities.  A^piv  is  the  number  of  PIV-sampled  points  on  the  centerplane. 

The  results  are  summarized  in  Table  1.  The  column  labeled  “(DC)  (%)”  contains  time-averaged  values 
of  the  jet  duty  cycle.  Values  in  the  rightmost  column  are  measures  of  control  effectiveness  or  efficiency  (the 
ratio  of  fluctuating  velocity  reduction  to  required  control  input): 

*  _  (^rms) Control  (^rms)No  Control  /q\ 

^  ^  (iXT)  ■ 

Higher  values  of  ^  correspond  to  more  efficient  controllers.  Two  of  the  simple  closed-loop  controllers  reduce 
velocity  fluctuations  more  than  the  representative  open-loop  case  (50%  duty  cycle).  Based  on  the  ratio  ^ 
in  the  rightmost  column  of  the  table,  all  of  the  closed-loop  controllers  operate  more  efficiently  than  the 
open-loop  system.  The  runs  with  the  band-pass  and  low-pass  filters  perform  most  efficiently,  and  data  from 
the  run  with  the  band-pass  filter  are  used  to  construct  an  advanced  controller. 


Run  Description 

(ttrms) 

(DC)  (%) 

No  Control 

9.65 

- 

- 

Open-loop  Control 

4.14 

50 

0.11 

Closed-loop  Control  (Band  Pass) 

4.67 

33 

0.15 

Closed-loop  Control  (Band  Stop) 

3.86 

44 

0.13 

Closed-loop  Control  (Low  Pass) 

3.47 

42 

0.15 

Closed-loop  Control  (Very  Low  Pass) 

5.14 

35 

0.13 

Table  1.  Summary  of  results  from  the  wind  tunnel  runs.  ^  is  defined  in  (3)  and  is  a  measure  of  controller 
efficiency. 


Figure  7  compares  results  from  the  band-pass-filtered  closed-loop  run,  an  open-loop  run,  and  the  baseline 
run  with  no  control  jets:  time  histories  of  performance  output  (rms  values  of  fluctuating  velocity  and  control 
input  (jet  momentum  coefficient).  From  the  top  right,  top  left,  and  bottom  right  panels  of  Figure  8,  the 
separation  region  is  moved  aft  of  and  below  the  aperture  by  the  initial  closed-loop  controller  as  it  was  by 
the  open-loop  controller. 


VI.  Controller  Design 

Next,  data  from  the  initial  closed-loop  run  are  used  to  construct  an  advanced  controller.  The  advanced 
controller  is  designed  in  Clear  Science  Corp.’s  software,  SMARTflow,  as  described  in  the  companion 
paper,  “Feedback  Flow  Control  for  a  Pitching  Turret  (Part  I).”  The  control  objective  in  the  initial  runs 
was  to  minimize  pressure  fluctuations  from  the  single  pressure  sensor,  but  the  true  objective  is  to  minimize 
velocity  fluctuations  (urms)  in  a  portion  of  the  aperture  field  of  view  close  to  the  aperture  surface.  The  actual 
field  of  view  is  three  dimensional,  but  the  PIV  window  intersects  this  volume  on  the  turret  centerplane.  PIV 
measurements  are  used  to  construct  the  estimators.  They  are  not  available  to  the  controller  in  real  time  but 
are  used  to  evaluate  controller  performance. 

The  single  pressure  sensor  used  in  the  initial  closed-loop  runs  was  chosen  because  it  correlated  well 
with  pitch  angle:  pressure  variance  increased  as  the  angle  (and  extent  of  flow  separation)  increased.  The 
correlation  seems  to  be  due  primarily  to  the  sensor’s  location  upstream  of  the  jet  slots,  which  will  change 
as  the  yaw  angle  changes.  A  more  robust  controller  will  utilize  multiple  sensors  by  relying  on  a  general 
pressure-velocity  correlation  contained  within  the  measurement-based  estimator  as  described  in  “Feedback 
Flow  Control  for  a  Pitching  Turret  (Part  I)” — rather  than  any  specific  relationship  between  surface  pressures 
and  velocity.  The  measurement-based  estimator  in  the  more  advanced  controller  utilizes  f8  sensors:  Sensor 
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Figure  7. 


Comparisons  of  the  baseline  (no  control),  open- loop  control,  and  initial  closed-loop  control  runs. 


Nos.  1,  2,  3,  5,  6,  8,  11,  12,  13,  14,  15,  16,  17,  18,  25,  26,  29  and  30,  identified  in  the  right  panel  of  Figure  1. 
These  sensors  maximize  the  correlation  between  pressure  and  velocity.  Pressure  sensors  are  located  on  the 
aperture  itself  for  the  purpose  of  investigating  the  flows  and  developing  preliminary  controllers;  however, 
these  sensors  would  interfere  with  operation  in  real  applications,  and  hnal  designs  would  exclude  them. 

The  first  step  in  constructing  estimators  is  to  compute  a  set  of  POD  basis  functions  using  PIV  snap¬ 
shots.  The  left  panel  of  Figure  9  contains  the  POD  energy  distributions  in  the  reduced-order  model 
(ROM)  built  with  snapshots  obtained  from  the  initial  closed-loop  run.  Four  POD  modes  contain  more 
than  80%  of  the  perturbation  energy,  and  this  number  of  modes  is  chosen  for  the  estimators.  The  right 
panel  of  Figure  9  contains  comparisons  of  time  histories  of  rms  velocity  fluctuations  from  the  PIV  data 
(labeled  as  “RAW”)  and  time  histories  obtained  by  projecting  the  PIV  snapshots  onto  the  four  POD  modes: 
Wrms(i)  ~  \/o'{n){t)a(^n){t)j  «  =  1, 4.  Even  with  only  four  modes  in  the  models,  the  approximations  based  on 
projections  are  very  accurate. 

The  next  step  is  to  use  the  POD  models  to  construct  the  measurement-based  and  dynamical  estimators 
described  in  “Feedback  Flow  Control  for  a  Pitching  Turret  (Part  I)” — both  of  them  linear.  Figure  10 
contains  time  histories  of  the  POD  coefficients  from  the  dynamical  estimators  (D-EST),  the  measurement- 
based  estimators  (M-EST),  along  with  the  values  obtained  by  projecting  the  PIV  snapshots  onto  the  four 
POD  modes.  From  the  top  left  panel  of  the  Figure  10,  predictions  from  both  types  of  estimator  are  accurate  in 
predicting  the  trends  in  the  first  POD  coefficient,  and  this  mode  contains  more  than  75%  of  the  perturbation 
energy.  Note  how  much  less  noise  there  is  in  the  dynamical  state  estimates  than  in  the  measurement-based 
state  estimates.  This  will  be  exploited  with  the  Kalman  filter. 

Again,  the  control  objective  in  the  wind  tunnel  is  to  minimize  velocity  fluctuations  in  the  aperture  field  of 
view  (or  at  least  over  the  area  of  intersection  between  the  PIV  window  and  field  of  view).  Figure  11  contains 
approximations  of  Urms  from  the  dynamical  and  measurement-based  estimators,  along  with  values  from  the 
projected  POD  coefficients.  From  Figure  11,  the  dynamical  estimator  is  able  to  accurately  track  the  velocity 


8  of  14 


American  Institute  of  Aeronautics  and  Astronautics 


X/D 


lu|/u  ■  I  I  I  I  I  I  I  I  I  I  I  rr^ 

n'innnon>incncn'Tnonn<in'i'i'inHO<i>i<ic<ic 


0.0  0.2  0.4  0.6  0.8 

X/D 


|u|/U„ 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0  1.1  1.2  1.3  1.4  1.5  1.6 


X/D 


|u|/U„ 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0  1.1  1.2  1.3  1.4  1.5  1.6 


Pressor  Sensor  #1 

0.2 1 - ^ ^ - 


0 


-0.2' - ^ ^ ^ ^ - 1 

5  10  15  20  25  30 

Pressor  Sensor  #14 


Pressor  Sensor  #28 


Q  2I _ ,  1  iirr  ^T..nprni-- III . .M '  71  T^ir'  M||  'I  f'l"  l[»l[r"  ”|r'l  Iff  1  "l  f!  111 

■  5  10  15  20  25  30 

Pressor  Sensor  #7 


. . . 

5  10  15  20  25  30 

Angle  of  Attaok 


Time  (sec) 


Figure  8.  Velocity  contours  and  streamlines  in  the  PIV  centerplane  from  the  run  with  dynamic  pitching  and 
simple  closed-loop  control  (Band-Pass)  at  dynamic  pitch  angles  of  115°  (increasing)  in  the  top  left  panel,  120° 
in  the  top  right  panel,  and  115°  (decreasing)  in  the  bottom  left  panel.  Time  histories  of  five  surface  pressures 
on  the  turret  centerline  from  the  dynamic  run  in  the  bottom  right  panel. 


Figure  9.  Energy  distribution  in  the  POD-based  ROM  (left).  Time  histories  of  Wrms  from  the  PIV  data  (RAW) 
and  approximations  obtained  by  projecting  the  PIV  snapshots  onto  the  first  four  modes  in  the  model  (PROJ) 
(right). 
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Figure  10.  Time  histories  of  the  first  four  POD  coefficients  from  the  dynamical  estimator  (D-EST),  the 
measurement-based  estimator  (M-EST),  and  valnes  obtained  by  projecting  the  PIV  snapshots  onto  the  four 
POD  modes  of  the  model  (PROJ). 


fluctuations.  The  measurement-based  estimator  exhibits  relatively  high  noise  levels — even  when  the  pressure 
signals  are  filtered — but  the  dynamical  estimator  and  Kalman  filter  described  in  “Feedback  Flow  Control  for 
a  Pitching  Turret  (Part  I)”  are  designed  to  address  this  issue. 

With  performance  output  defined  as  spatially  integrated  values  of  Urms(t),  the  controller  will  utilize 
approximations  from  the  estimators,  and  the  PI  regulator  described  in  “Feedback  Flow  Control  for  a  Pitching 
Turret  (Part  I)”  becomes 


=  -  0)  +  k2  (^£n(T)x„(T)  -  0)dT 

Jo 


(4) 


where  u  is  the  suction  valve  duty  cycle  percentage  and  the  objective  is  to  reduce  fluctuating  velocity  to  zero. 
The  gains  (^1,2)  are  determined  heuristically,  which  can  be  a  time-consuming  process  involving  a  number  of 
wind  tunnel  runs.  Alternatively,  a  Linear  Quadratic  Regulator  (LQR)  automatically  determines  gains  that 
minimize  the  following  cost  function: 


1  F 

{XnXn  +  ru^)dt. 

J  Jo 


(5) 


The  parameter  r  weights  the  penalty  assigned  to  the  control  input  (u). 

Next,  the  dynamical  estimator  from  “Feedback  Flow  Control  for  a  Pitching  Turret  (Part  I)”  is  written 
as 

X  —  Gq  +  Aix  +  bou  +  b2U  +  Cw  and  w  =  Hw  (6) 
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TIME 


Figure  11.  Approximations  of  Urms  from  the  dynamical  estimators  (D-EST),  measurement  estimators  (M-EST), 
along  with  values  obtained  from  projected  POD  coefficients  (PROJ). 


where 


C=ib3  64),  w  =  and  H  =  (7) 

uj  is  the  pitching  frequency,  defined  in  (1).  Defining  an  augmented  state  estimate  a.s  z  =  {x,u,w)'^  and 
minimizing  the  cost  function  in  (5)  produces  the  following  control  law: 

il  =  (8) 


The  dynamical  estimator  without  the  Kalman  filter  were  imported  from  SMARTflow  into  Matlab,  and 
the  following  gains  were  computed: 

j^lqr  ^  (•_o.2746  0.5013  0.0482  -  0.0949  0.0064  -  1.3309  -  3.2540)^  x  10^.  (9) 


The  gains  in  (9)  were  computed  after  reducing  the  weighting  factor  r  in  (5)  to  its  lowest  possible  value, 
which  effectively  imposes  no  penalty  on  the  control  input. 


VII.  Closed-loop  Control  Run  with  an  Advanced  Controller 

Three  compensators  (LQR  regulators  with  Kalman  filters)  were  designed  using  three  different  values  of 
the  process  covariance.  As  the  process  covariance  decreases,  the  controller  relies  more  on  the  dynamical 
estimator  and  less  on  the  measurement-based  estimator — see  “Feedback  Control  for  a  Pitching  Turret  (Part 
I)”  for  a  definition  of  the  term.  From  Table  2,  the  process  covariance  is  highest  in  Run  1  and  lowest  in  Run 
3.  Measurement  error  covariance  is  3.3  x  10“"^  in  all  three  runs. 

Figure  12  contains  time  histories  of  fluctuating  velocity  (middle  panel)  and  jet  momentum  coefficient 
(bottom  panel)  from  the  three  control  runs,  along  with  baseline  histories  from  the  uncontrolled  run.  Com¬ 
paring  the  blue,  pink,  and  red  curves  in  the  middle  panel,  it  is  clear  that  heavy  reliance  on  the  dynamical 
estimator  results  in  poorer  performance.  The  best  control  is  achieved  in  Run  1  where  the  measurement-based 
estimator  dominates.  From  the  pink  and  red  curves  in  the  bottom  panel  of  Figure  12,  the  Kalman  filter 
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Quantity 

Symbol 

Run  1 

Run  2 

Run  3 

Signal-to-Noise  Ratio 

SNR 

5.0 

5.0 

5.0 

Process  covariance 

Wp 

3.3  X  10-1 

3.3  X  10-1 

3.3  X  10-1' 

Measurement  covariance 

Vp 

3.3  X  lO-'i 

3.3  X  10-1 

3.3  X  10-1 

Table  2.  Parameter  settings  for  the  wind  tunnel  control  runs. 


reduces  noise  in  the  control  input  significantly  when  the  dynamical  estimator  is  heavily  weighted;  however, 
the  blue  curve  indicates  that  substantial  filtering  is  effected  even  in  Run  1  when  the  dynamical  estimator 
is  weighted  much  less:  compare  the  relatively  smooth  blue  curve  in  the  bottom  panel  of  Figure  12  with  the 
noisy  pink  curve  in  the  bottom  panel  of  Figure  7. 

From  the  bottom  panel  of  Figure  12,  the  control  input  modulates  with  pitch  angle  in  all  three  runs  with 
the  advanced  controllers.  The  phase  shift  between  pitch  angle  and  input  varies  between  the  three  runs,  but 
the  turret  motion  is  clearly  driving  the  input,  and  modulation  in  the  jet  amplitude  is  not  small.  Conversely, 
the  pink  curve  in  the  bottom  panel  of  Figure  7  indicates  little  or  no  correlation  between  control  input  from 
the  simpler  controller  and  pitch  angle:  modulations  in  the  jet  momentum  coefficient  are  essentially  random. 


Figure  12.  Comparison  of  pitching  runs  with  and  withont  control:  results  from  the  new  controller. 


Table  3  summarizes  results  from  the  various  runs:  time-averaged  values  of  fluctuating  velocity,  control 
input  {{DC)),  and  controller  efficiency  (^) — defined  in  (3).  Higher  values  of  ^  correspond  to  more  effi¬ 
cient  controllers.  Again,  the  run  with  a  steady  suction  jet  operating  at  50%  duty  cycle  is  chosen  as  the 
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representative  open- loop  case.  It  effects  the  largest  reduction  in  fluctuating  velocity  (57%);  however,  it  is 
significantly  less  efficient  than  the  two  most  efficient  closed-loop  controllers.  The  most  efficient  simple  con¬ 
troller  (band-pass)  reduces  fluctuating  velocity  by  52%  and  the  most  efficient  advanced  controller  (Run  1) 
by  51%.  Comparing  the  bottom  left  panels  of  Figures  6,  8,  and  13,  separated  flow  is,  on  average,  closer 
to  the  aperture  at  a  pitch  angle  of  115°  (decreasing)  in  the  advanced-controller  run  (Run  1)  than  it  is  in 
the  open-loop  and  simple  closed-loop  control  runs — reflecting  the  slightly  higher  time-averaged  fluctuating 
velocity  in  the  advanced-controller  run.  Comparing  the  bottom  right  panels  of  Figures  6,  8,  and  13,  the 
behavior  of  fluctuating  pressures  from  the  hve  sensors  on  the  turret  centerline  is  qualitatively  similar  in  the 
open-loop,  simple  closed-loop,  and  advanced  closed-loop  runs.  From  the  values  of  ^  in  Table  3,  the  advanced 
controller  in  Run  1  is  13%  more  efficient  than  the  most  efficient  simple  controller  and  55%  more  efficient 
than  open-loop  control. 


Run  Description 

(Wi-ms) 

{DC)  (%) 

No  Control 

9.65 

- 

- 

Open-loop  Control 

4.14 

50 

0.11 

Simple  Closed-loop  Control  (Band  Pass) 

4.67 

33 

0.15 

Advanced  Closed- loop  Control  (Run  1,  Wp 

>Vp) 

4.69 

30 

0.17 

Advanced  Closed- loop  Control  (Run  2,  Wp 

=  Vp) 

5.64 

40 

0.10 

Advanced  Closed-loop  Control  (Run  3,  Wp 

<Vp) 

7.97 

41 

0.04 

Table  3.  Summary  of  results  from  the  wind  tunnel  runs  with  open-  and  closed-loop  control.  ^  in  the  rightmost 
column  is  defined  in  (3)  and  is  a  measure  of  control  efficiency — the  ratio  of  fluctuating  velocity  reduction  to 
required  suction. 


VIII.  Conclusion 

Closed-loop  systems  for  controlling  flow  separation  and  turbulence  above  a  pitching  turret  have  been 
developed  and  tested  through  a  series  of  control-in-the-loop  CFD  simulations  and  wind  tunnel  runs.  The 
wind  tunnel  tests  have  been  described  here;  controller  designs  and  the  computational  tests  are  described 
in  the  companion  paper,  “Feedback  Flow  Control  for  a  Pitching  Turret  (Part  I).”  Levels  of  flow  separation 
and  turbulence  change  as  the  turret  pitches,  and  the  systems  successfully  reduce  these  levels  in  the  dynamic 
environment  through  feedback  control — meeting  the  objective  of  control  in  the  presence  of  disturbances  over 
a  range  of  operating  conditions.  Another  objective  is  the  minimization  of  control  input:  achieving  threshold 
levels  of  control  with  less  actuator  energy.  Closed-loop  systems  in  the  wind  tunnel  runs  are  as  much  as  55% 
more  efficient  than  a  representative  open-loop  system  as  measured  by  the  ratio  of  reduction  in  fluctuating 
velocity  to  required  control  input. 
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Figure  13.  Velocity  contours  and  streamlines  in  the  PIV  centerplane  from  the  run  with  dynamic  pitching  and 
closed-loop  control  (Run  1)  at  dynamic  pitch  angles  of  115°  (increasing)  in  the  top  left  panel,  120°  in  the  top 
right  panel,  and  115°  (decreasing)  in  the  bottom  left  panel.  Time  histories  of  five  surface  pressures  on  the 
turret  centerline  from  the  dynamic  run  in  the  bottom  right  panel. 
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